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We discuss various aspects of most general multisupport solutions to matrix models in the 
presence of hard walls, i.e., in the case where the eigenvalue support is confined to subdomains of 
the real axis. The structure of the solution at the leading order is described by semiclassical, or 
generalized Whitham-Krichever hierarchies as in the unrestricted case. Derivatives of tau- functions 
for these solutions are associated with families of Riemann surfaces (with possible double points) 
and satisfy the Witten-Dijkgraaf-Verlinde-Verlinde equations. We then develop the diagrammatic 
technique for finding free energy of this model in all orders of the 't Hooft expansion in the reciprocal 
matrix size generalizing the Feynman diagrammatic technique for the Hermitian one-matrix model 
due to Eynard. 

The matrix models and their so-called multisupport (multicut) solutions became important again 
after the studies in M = 1 SUSY gauge theories P, [2| and due to Dijkgraaf and Vafa [3], who proposed 
to calculate the low energy superpotentials using the partition function of multicut solutions. These 
solutions, known since early 90th (see, e.g., [HE]), were revived by Bonnet, David and Eynard 

Dijkgraaf and Vafa proposed to consider the nonperturbative superpotentials of N = 1 SUSY gauge 
theories in four dimensions arising from the partition functions of the one-matrix model (1MM) in the 
leading order in 1/N, N being the matrix size. The leading order (of the 't Hooft 1/iV-expansion) of 
the matrix model is described by the semiclassical tau-function of the so-called universal Whitham- 
Krichever (WK) hierarchy [7] (see also [H1IH1E3; the details about one-matrix and two-matrix cases 
see in ^2] and JS])- At the same time, matrix integrals beyond large- N limit are tau- functions of the 
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hierarchies of integrable equations of the KP/Toda type [Tl]. For the planar single-cut solutions, the 
matrix-model partition functions become r-functions of the dispersionless Toda hierarchy, one of the 
simplest example of the WK hierarchy. 

One may also consider more general solutions to matrix models, identifying them with generic 
solutions to the loop (Schwinger-Dyson, or Virasoro) equations ^3], or the Ward identities satisfied 
by matrix integrals ^H]. We discuss an interesting class of multi-cut, or multi-support, solutions to 
the loop equations possessing multi-matrix integral representations |3 El EEEI E3- These solutions are 
associated with families of Riemann surfaces and form a sort of a basis in the space of all solutions to the 
loop equations |181ll9j (like the finite- gap solutions form a similar basis in the space of all solutions to an 
integrable hierarchy). They can be distinguished by their "isomonodromic" properties — switching on 
higher matrix model couplings or 1 /^-corrections does not change the family of Riemann surfaces, but 
just reparameterizes the moduli as functions of these couplings. This property is directly related to that 
the partition functions of these solutions are semiclassical tau-functions (also called the prepotentials 
of the corresponding Seiberg-Witten-like systems). 

The geometrical properties of these solutions were extensively studied in [201 121j . and the com- 
prehensive procedure for calculating 1 /^-expansion terms was developed in [23]. In this paper, 
we give the review of methods of these papers and extend them to the case of matrix models with 
hard walls (hard edges). Those are models in which we introduce rigid restrictions on the domain of 
admissible eigenvalues of the potential. It is well-known that the behavior of the eigenvalue distri- 
bution changes dramatically in the presence of such walls: it has inverse square-root singularity as 
approaching the wall (see, e.g., Szeg0 |2lj). 

Our aim in this paper is to show how the solutions with hard walls can be nicely incorporated in the 
general setting of solutions to 1MM. Actually we show that all the basic ingredients of 1MM solutions, 
that is, the WK hierarchy, the WDVV equations, and the diagrammatic technique for evaluating 
corrections in 1/N, turn out to be common for solutions with and without hard walls and only the 
hyperelliptic Riemann surface representing the spectral curve is changed. Technical solutions however 
apparently develop distinctions, which we indicate in what follows. 

Recently, models with hard walls in the planar limit of two-matrix models were constructed |251l26j . 

For the generality of presentation, we use the term 1MM to denote a Hermitian one-matrix model 
possibly with hard walls. 

In Sec. n we describe the general properties of multi-cut solutions of matrix models. We derive 
the basic set of constraints on the matrix model partition functions and present them in the form of 
the master loop equation. 

Before presenting the general procedure of solving the master loop equation order by order to find 
the loop mean (one-point resolvents) of the matrix model potential, we present the special investigation 
of just the leading (planar) approximation. We prove that the free energy of the 1MM in the planar 
(large- N) limit coincides with the prepotential of a Seiberg-Witten-Whitham-Krichever theory. 



2 



The corresponding WK hierarchy is formulated in terms of Abelian differentials on (a family 
of) Riemann surfaces. This implies that the main quantities in matrix models can be expressed 
in geometrical terms and allow calculating derivatives of the matrix model free energy. Indeed, we 
demonstrate in Sec. [21 that the second derivatives of the 1MM free energy can be expressed through 
the so-called Bergmann bi-differential. On this language we derive the general relation satisfied by the 
matrix model potentials when differentiating w.r.t. Whitham times. We then come to third derivatives. 
Indeed, it is typical (but not compulsory) for WK tau-functions to satisfy a set of (generalized) 
Witten-Dijkgraaf-Verlinde-Verlinde (WDVV) equations [23 [2E1 [2U that are differential equations 
with respect to Whitham times involving third derivatives. These equations are usually considered an 
evidence for an underlying topological theory. In Sec. [21 we prove that the semiclassical tau-function 
of the multi-support solutions to matrix models satisfies the WDVV equations in the case of general 
1MM solution, i.e., in the case of arbitrary number of nonzero canonical times, which include the 
times of the potential and the occupation numbers (filling fractions), which indicate the portions of 
eigenvalues that dwell on the related intervals of the eigenvalue support. The hard- wall parameters, 
although being out of the set of canonical variables, affect nevertheless the dimension of the set of 
nondegenerate WDVV equations. This completes an interpretation of the results of ^ [21 in terms of 
semiclassical hierarchies. The WDVV equations follow from the residue formula and associativity of 
some algebra (e.g., of the holomorphic differentials on the Riemann surface) |29| I3()|. We present the 
proof of the residue formula and extra conditions following 20 . Analogous results hold true in the 
two-matrix model 3 1 

In Sec. [31 we describe the diagrammatic procedure that allows evaluating order by order all the 
corrections in 1/N in the 1MM with hard walls. In the same section we describe the procedure of j2B] 
to integrate the obtained answer to obtain the expression for the free energy itself. This procedure 
works in all genera except the subleading (torus) approximation, which, correspondingly, deserves the 
special investigation (as in the case of the standard 1MM). We therefore devote the special section @] 
to the calculation of the free-energy subleading correction. Note again that analogous formulas were 
obtained in the two-matrix model case by Eynard, Kokotov, and Korotkin 32 . 



1 Matrix models and Riemann surfaces 

1.1 Matrix integrals and resolvents 

Consider the 1MM integral 

JNxN;E 

where V(X) = J2k>itkX k , h = t8> is a formal expansion parameter, the integration goes over the 
N x N matrices, DX oc Yli<j d^tXydlmXij J\ i dMXa, and we also assume that the integration domain 
E for the eigenvalues of the matrix X comprises a number of (possibly half-infinite) intervals, E = 
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u|=i[°2/?-i) 02/3]- The topological expansion of the Feynman diagrams series is then equivalent to the 
expansion in even powers of h for 

00 

T = T{h,to,t 1 ,t 2 ,...) = J2 n2h:F h, (2) 

h=0 

Customarily to = fi>N is the scaled number of eigenvalues. We assume the potential V(A) to be a 
polynomial of the fixed degree m + 1, with the fixed constant "highest" time t m+ \ = 1. 

The averages, corresponding to the partition function Q are defined as usual: 

</(*)> = \\ DX f(X) exp f-itr V(X)) (3) 
£ JnxN;e V n J 



and it is convenient to use their generating functionals: the one-point resolvent 

;trx* 



k=0 



as well as the s-point resolvents (s > 2) 

W(h,....K)-V t ^ • ^11°" = (tr ^ ■ ■ ■ tr ^ ) (5) 

fel ... /b s =l • • • A s \ A l A A s -A- I conn 

where the subscript "conn" pertains to the connected part. 

These resolvents are obtained from the free energy T through the action 



dV(X s ) dV(X s ^) dV(X 
of the loop insertion operator 



dv(x s ) dv(x 8 -!) av(\i) 

d d d 

W(\i), (6) 



a yJ-± ( 7) 



9V(X)- fyX^dty 

Therefore, if one knows exactly the one-point resolvent for arbitrary potential, all multi-point resolvents 
can be calculated by induction. In the above normalization, the genus expansion has the form 

00 

W(X 1 ,...,X s ) = Y / K 2h W h (X 1 ,...,X s ), s>l, (8) 

h=0 

which is analogous to genus expansion (J2J). 
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1.2 Master loop equation 

We begin with deriving conditions on integral (^Q) coming from its invariance w.r.t. changing the 
integration variables. For this, we first to the eigenvalue representation assuming 

eV* = / DXe-%* v m ^ TT ^A 2 (A)e-^^>^) (9) 

JE({a p }) JE{{a p }) H 



up to the volume of the unitary group U(to/fb). Here, as usual, we let A (A) denote the Vandermonde 
determinant of variables Aj. In the latter expression, we perform the change of variables 

A 4 = A 4 + eA| +1 , P >-1, (10) 

obtaining 

to /ft 

e^= TT d(A l + ,Ar 1 )A 2 (A + e A^ 1 )e-^t / >(A l+£ Ar 1 ). 



i=l 



Segregating the coefficient by e in a standard way and equating it to zero, we obtain the Virasoro 
conditions to be satisfied by integral Q (see [35] ) 

^ H ^ otk otp-k £? x ot k+p I 



where by definition 

— Z = - — Z or— J 7 = -t 



In terms of J-, Eqs. (fllj) become 

^ p da/j ^ dt k dt p - k ^ Q dt k dt p _ k ^ dt k+p 

and contracting these Virasoro conditions with A~ p ~ 2 , we obtain the master loop equations ^Hj of the 
1MM with hard walls Ql]: 

tl J —^ +< ( ^^W(0 = W(\) 2 + h 2 W(\,\). (13) 
^ A - a/? da/? JC B 27r« A — £ 

We introduce the linear integral operator 

Kf{\) = f = [V(A)/(A)]_ (14) 
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projects onto the negative powers 1 of A. Hereafter, Cx> is a contour encircling all singular points of 
W(£), but not the point £ = p. Using Eq. @, one can express the second term in the r.h.s. of loop 
equation (|T3|) through W(\), and Eq. (|T3)) becomes then an equation for one-point resolvent Q. 

Substituting genus expansion (jHJ) in Eq. (jl.Sj) . one finds that Wh(X) for h > 1 satisfy the equation 



E I TT^ + [ K ~ 2W °( X )) W hW = E W h ,{\)W h - h ,{\) + ——W h ^(X), (15) 

^X-apdap \ ' ^ 0V\A) 

In Eq. (|T5*|) . Wh(A) is expressed through only W^iX) for which hi < h. This fact allows developing 
the iterative procedure. 

Note that the loop equation contains now explicit dependence on the points ap. Although the 
term proportional to the derivative of Th seems not to be universal, in what follows we find that these 
terms are important to keep the picture self-consistent. 

We use the asymptotic condition (which follows from the definition of the matrix integral) 

Wfc(A)U-oo = j5 hfi + 0(1/A 2 ), (16) 

and manifestly solve (|13|) for genus zero. Then, one could iteratively find Wh(X) thus restoring the 
corresponding contributions in the free energy by integration, since 

W *W = -W(xf*> h>l. (17) 

1.3 Solution in genus zero 

In genus zero, loop equation (fT3)) becomes 

E + / ££^(0 = [V'(X)W (X)]_ + =*S=m = mX )f, (18) 

where the polynomial P 9 _i(A) is such that 

y _^^£ = ^-i( A ) ( i9) 

^X-ap dap 11/3=1 ( A - a p) 



1 In order to prove it, one suffices to deform the integration contour to infinity to obtain 

= V'(\)f(\) - [V'(\)f(\)] + = [V'(X)f(X)]_ 
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In what follows, for simplicity, we let P r (\) denote any polynomial of A of order r (which is 
indicated by the subscript). For instance, we can write P r (X) + P U (X) = P W (X), where w = max(r, u). 

In order to solve Eq. (|18|). note that 

[V'(X)W (X)]_ = V'{X)W Q {X) - [V'(X)W (X)] + . (20) 
Due to ()16|) . the last term in the r.h.s. is a polynomial of degree m— 1, m being the degree of V'(A), 

P m -i(A) = -[V(A)Wb(A)] + = -jf ^X^WbCO- (21) 



That is, we have 



w*(x)-v'(x)w (x) = p m „ 1 (x)+ P(? - l(A) p ^-^) 



xiUxix-ap) n^=i(A-^)- 

Then, the solution to (|18|) is 



1, l^(APn^(A- gri+ 4P m+1 - lW 

2 2 'nLi( A -°-3) 



where the minus sign is chosen in order to fulfill asymptotics (|16|) . We also assume that the general 
solution in the large- iV limit depends only on t among q hard wall points an. Then, from (|19|) . we 
have that the corresponding pole disappear and does not enter expression ([22)1 . We therefore present 
the one-point resolvent in the form 

Wo{X) = \{V'{X)-y{X)), (23) 
therefore introducing a new function y{X) satisfying the equation 

y(Xf = V\Xf + A j t m+t ' l{X \ ( 24 ) 

1 1/3=1 ( A - a /3) 

and determined on a hyperelliptic Riemann curve. For the generic potential V(A) with m — * oo, this 
curve may have an infinite genus, but we are going to consider solutions with only finite number n of 
cuts. Endpoints of these cuts are t points ag (t > q), which we enumerate by /3 = 1, . . . , t, and s new 
endpoints fj, a , a = 1, . . . , s. Apparently, s + t = In. We then present y(X) in the form 



,2,^ _ U s a =i(x-fi a ) 
UB=i(X-ap) 



y(X) = M(X)y(X), and f{X) = ^ (25) 
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with all fi a , ap distinct. In what follows, we assume M(X) to be a polynomial of degree m — 2jp, 
keeping in mind that both s and t are always finite and fixed, while m > n can be chosen arbitrarily 
large. By convention, we set y( A )U-+oo ~ A^ s_ ^' 2 , and due to (|23j). 

y'(A)=y(A) + ^ + 0(l/A 2 ) (26) 

at large A. 

In comparison with the standard 1MM situation, we may encounter an ambiguity. Indeed, as y(A) 
is proportional to A*- S_i ^ 2 as A — ► oo, multiplying Wo(£) by ^_->p and integrating at infinity may 
produce a nonzero result for t > s. For M(A), we then have 



s — t 
2 



m(a)=/ |i ^ : ' ^ n ( a -^)- ( 27 ) 

Here Pt-s (6) is a polynomial and 6j are zeros of M(A). In our consideration, we systematically 

2 

avoid any explicit dependence on 6j. Instead, the relevant objects on which the free energy depends 
are moments of the general matrix model potential (see |M4| I35[ I17j ) representing (p — l)st 

derivatives of this polynomial at the branching points that are endpoints of intervals of eigenvalue 
distribution, 



m (p) = ( _f!___M(A) 

(p-l)\ \d\P- 1 v ; 



, e.g., MP = M(m). (28) 

Among endpoints m of these intervals there are hard wall points ag, (3 = l,...,t (t < q) and a 
number of new "dynamical" endpoints fj, a , a = 1, . . . , s with s + 1 = In. Often (but not always) these 
points can be considered on the equal footing, as in formula (|28|). and we reserve the unified notation 
Hi with the Latin subscript i = 1, . . . , 2n to indicate any branching point, whatever origin it has. 

Directly from definition (|28|) follows the representation for the moments 

M? ] =i £ ' ' 2 , P>1. (29) 



Inserting solution (|25|). (|27|) into (|23|) and deforming the contour, we obtain the planar one-point 
resolvent with an n-cut structure, 

The contour Cd of integration here encircles the finite number n of disjoint intervals 

n 

V = [J\P>2i-l, fl>2i\, Ml < V2 < ■■■ < M2n- (31) 

i=l 
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We now address the question how many free parameters do we have. In (|3Uj) we have m times t& 
of the potential V^(£) (recall that the highest time is the unity), ^ — 1 coefficients of the polynomial 
Pt-s (£) (the constant term does not contribute) and s new endpoints \x a . However, we must take into 

2 

account asymptotic conditions which impose the restrictions 

- 2f <L=t -tl) , 10 = 0,1,...,——, (32) 



if s > t, or 



2 



d£ "— _ . 



2 



7 2to, where Pt-s (A) = at- s A 2 +..., (33) 

ICd 2m y(£) 2 2 

for s < t. If we do not count the hard wall parameters, in the first case, we have m times of the 

potential, s points /x a , one zero time to minus + 1 restrictions, that is m+ = m + n parameters. 

In the second case, we have m times of the potential, s points /J, a , one zero time to, ^ relevant 

coefficients of the polynomial Pt-s (A) and one restriction, which gives again m + s ^ = m + n free 

parameters. That means that, except the time to, we have an n — 1-dimensional space of parameters. 

We can arbitrary choose coordinates on this space, but it turns out that there exists a distinguished 

set of n — 1 independent variables that parameterize solutions to the loop equations [HI El , 

5? (34> 

where Ai, i = 1, ... ,n — 1 is the basis of ^4-cycles on the reduced hyperelliptic Riemann surface of 
genus h = n — 1 

s t 2n 

f = n ( a - n ( a - = ik a - ( 35 ) 

«=1 /3=1 i=l 

(We may conveniently choose the ^4-cycles to be the first n — 1 cuts, see Fig.l.) Note that we do 
not distinguish the branching points fj,i in the above expression by their origin: both "dynamical" 
branching points fi a and hard walls ap enter y(X) on the equal footing, in contrast to the standard 
1MM case, y(A) 7^ y(A). However, in what follows, y(\) enters some of the constructions of the paper. 

Given the basis of A-cycles, we also choose the conjugated basis of B-cycles with the intersection 
form AioBj = 6™. Besides canonically conjugated A- and -B-cycles, we also use the linear combination 
of -B-cycles: P>i = Bi — Bi+i, P> n -\ = Pn-\- Therefore, B-cycles encircle the nearest ends of two 
neighbor cuts, while all -B-cycles goes from a given right end of the cut to the last, nth cut. For the 
sake of definiteness, we order all points \ii in accordance with their index so that ^ is to the right of 
Hj if i > j. 

1.4 Matrix eigenvalue picture 

The variables Si find a "physical" interpretation in the semi-classical picture of matrix eigenvalues 
(Coulomb gas) in the limit where their number (and, therefore, size of the matrix) goes to infinity. 
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Si s*2 53 Sg to — J2 Si 

A, B i A 2 B 2 As Ag B,=B, An 



Figure 1: Structure of cuts and contours for the reduced Riemann surface. 

Let us first introduce the averaged eigenvalue distribution 

f N 1 

PW = Jj E W A - A *)> = J™ ( A - * e ) - W ( X + ^)) ( 36 ) 

i 

with Aj being eigenvalues of the matrix X. In the planar limit, this quantity becomes 

po(X) = ^-lim(w (X-ie)-W (X + ie)) = ^lmy(X) (37) 
and it satisfies the equation 2 

/ E^L dX = l V '{x), VpeV. (38) 
Jv x — X 2 

This averaged eigenvalue distribution becomes the distribution of eigenvalues in the limit where their 
number goes to infinity. Then, matrix integral Q becomes 

eM = | E n^^^ (/y(A)e(A) " /e(A)e(A ' )l0g|A_A ' |dAdA ' ) = J E Y[d Xl e^ (39) 
where we introduced the eigenvalue distribution 

<?(*) = ^ ( 4 °) 

i 

Now, in the limit of large N, one can use the saddle point approximation to obtain the equation for 
q(X). Since the variable to plays the role of the (normalized) total number of eigenvalues, 

If 1 

to = 4vri fc y ^ dX = ~ 2 res °°( ydA ) ( 41 ) 



2 Indeed, by definition 

W (X) 



, dX = 0. 

x — A 

Now, using ll2Mt and definition l|36^ and puffing out the contour from infinity, one easily comes to this equation. 
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and the support V of g consists of n segments T>i, we impose the constraint 

g(X)d\ = t , (42) 



and, following jSJ Ej, fix the occupation numbers of eigenvalues in each of the segments, Si (|34JI . 
i = 1, . . . ,n — 1. We assume the occupation number for the last, nth cut to be to — J2?=i Si = S n . 3 
(Obviously, no new parameters Si arise in the one-cut case.) We formally attain this by introducing 
the corresponding chemical potentials (Lagrange multipliers) Ho and IT, i = 1, . . . , n — 1, in the 
variational problem for the free energy, which therefore becomes in the planar limit 



S eS [g; Si,t , t k ] = - V{X)g(X)dX + / / g(X) log |A - A'| g(X')dXdX' + 

Jv J Jv 



-n g(X)dX - *o) - E Hj g(X)dX - S 



(43) 



while the saddle-point equation becomes 

2 J g(X) log |x - X\dX = V(x) + 11$ + n , Vx G A, (44) 

and its derivative still coincides with (|38|). 

Therefore, with generic values of the constants LTj, g c (A) gives the general solution to (|3*H)l (or the 
planar limit of the loop equation) : these constants describe the freedom arising when solving the loop 
equation. However, in the matrix model integral (where there are no any chemical potentials) one 
would further vary To w.r.t. LTj to find the "true" minimum of the eigenvalue configuration, 

f -0, V, (45) 

This is a set of equation that fixes concrete values of Si and Ilj in the matrix integral. 
Let us now calculate the derivative of Fo (j4*3|) w.r.t. Si. From (|43j) . one has 

dS eS 



dSi 



e=po 



v 



The expression in the brackets on the r.h.s. of (|46[) is almost a variation of (|43[) w.r.t. the eigenvalue 
density, which is 



q ^ S'eff 



8p(X) 



e=po 



V(X) -2 [ dX' log(A - A')po(A') + Hi + n c 



for XGDiCV. (47) 



3 It is sometimes convenient to consider S n instead of to as a canonical variable. However, in all instants we use S n 
we specially indicate it for not confusing S n with the "genuine" filling fraction variables Si, i = 1, . . . , n — 1. 
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It is therefore a step function h(X), which is constant equal to d = — ITg — H", on each cut Aj. One 
then has 



dF dS e{i [g] 



dS 4 dS, 



Q=P0 J u ""I ' j = l 



OS 

j'=i 

In particular, 



E0gg7 = -<. + C„ = ni. (48) 



In [3] it was proved that the difference of values of Ilj on two neighbor cuts is equal to 4 

Ci+i - Ci = 2 / Wo(A)dA, (50) 



i.e., 



Ili = (Ci+l - Ci) + (Ci+2 - Ci+l) + • • • + (Cn-l - Cn-2) + (Cn - Cn-l) 



y(A)o!A = * y(A)dA. (51) 

B,UB i+ iU...UB 9 JB, 

The planar limit of the free energy can be obtained by substituting the saddle point solution g 
into (1331) : 

If 11 nl 

•To = S eS [q c ] = -- / V(X)g c (X)dX + -n t + ViSi. (52) 
2J V 2 2 , =1 

In the paper, we choose the solution to the loop equation with fixed occupation numbers, (|34[). 
Note that fixing the chemical potentials (|48 [) -(|51 [) instead of (]45[). we obtain the solution with the 
interchanged A- and -B-cycles on the Riemann surface (|24|). However, .Fq is not modular invariant. 
Under the change of homology basis, J-q transforms in accordance with the duality transformations, 
see [311! • The higher-genus corrections become also basis-dependent: choosing Si or IL as independent 
variables, one obtains different expressions, say, for the genus-one free energy, see Sec. 4.3. 

Note that the presence of hard walls has had no effect on consideration in this section because 
it only changes asymptotic behavior of the function y(A) near the corresponding branching points: 
while for "dynamical" branching points fj, a the eigenvalue distribution undergoes the semicircular 
Wigner's law, it is known (see, e.g., [21]) that in the vicinity of hard wall we observe accumulation of 

4 The simplest way to prove it is to define function h(X) outside the cuts: h(X) — V(X) — 2 J dA'log(A — X')po(X') 
and note that h'(X)\ x ^ v = 2W (X). 
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eigenvalues (zeros of the corresponding orthogonal polynomials) with the reciprocal square-root law. 
This perfectly matches with our answer (I25[) for y(A). 

In the next subsection, we demonstrate that the planar loop equation solution with fixed occupation 
numbers corresponds to a Seiberg-Witten-Whitham-Krichever system, just as in the case of ordinary 
1MM [20]. 



1.5 Free energy as the prepotential of SWWK system 

We now turn to studying geometrical properties hidden in the matrix- model solutions. Namely, we 
associate a Seiberg-Witten-Whitham-Krichever (SWWK) system with the planar limit of the matrix- 
model free energy. 



Matrix integral as a Seiberg Witten system. The family of Riemann surfaces is now the family 
of g = n — 1 reduced Riemann surfaces described by (|24|) or (|25[). This means that these Riemann 
surfaces contain no information about the additional polynomial M(A), which is present, however, 
through (|2l)|) in the differential dS. The role of SW differential is played by 

dS = y{X)d\. (53) 

We now consider its variation w.r.t. Si- Variations over moduli of the surface do not change the 
genus of the reduced Riemann surface as well as the highest degree of the additional polynomial 
M m _ n (A). We also consider both the times of the potential V^(A) and the hard wall parameters ap to 
be independent on the parameters Si and to, that is, we assume SV 1 /5Si = and consider the general 
variation 5, which varies the potential and the branching points /j, a (but not the hard wall parameters 
a/3 and the special variation 5s that leaves invariant the potential and the hard wall parameters and 
change only the moduli parameters Si. 

Using ([24)1 . (|25|). one obtains for the general variation 5dS: 5 



-dX, (54) 



SdS = 5 (M m _w(A)y(A)) dX = 

nU(A - n*)8M m ^(X) + iM m _^(A)<5ffi=i(A - ^ 



m 

where the polynomial expression in the numerator is of maximum degree m + ^- — 1 = m + n — 1 (since 
the highest term of M m s -t is fixed) and we have in the denominator the hyperelliptic coordinate 
(|35[). On the other hand, under 5s which does not alter the potential, we obtain from (|24[). (|25|) that 

5 s ±P m +t-i{X) 



5sds = -~- — ,r .. — Tdx 

1/3= 



2M m _ £ ^(A)y(A)n^=i(A-a /3 ) 



5 Note that the variation S differs nevertheless from loop insertion @ because the former does not change, by definition, 
the degree of the polynomial M(A). 
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- ^ gPm+ Mwl^ . (55) 
M m _^(A)y(A) V ; 



Because this variation is just a particular case of (|54|). we obtain that zeros of M a -t (A) in the de- 
nominator of l|55|) must be exactly cancelledby zeros of the polynomial <5s_P m +t_i(A) in the numerator, 
so the maximum degree of the polynomial in the numerator is — 2 = n — 2 (because, again, the 
highest-order term of P m+ t-i(X) is fixed by asymptotic condition (|26|) and is not altered by variations 
6s)- We then come to the crucial observation that the variation 5sdS is holomorphic on curve ()35|) . 
as it should be for the SW differential. 

Moreover, considering specific variation w.r.t. Si with condition (|34j) taken into account, we find 
that 

5 itj = S SiSj = S Si £ ^y(X) = £ ^«ay(A) 

and hence 

ddS Hj{\)d\ . 
-WK- = ~ m = (56) 
d^i y(X) 

with dwj being the canonically normalized fg^j du>i = holomorphic 1-differentials on the re- 
duced Riemann surface y(A). Here Hi(\) are polynomials of degrees n — 2 (the degree is exact because 
each such polynomial has exactly one zero on each cycle Aj for j ^ i). 

Rewriting now (|BT|) as = § B . dS, together with l|56|). we reconstruct the SW system [37| whose 
prepotential is the planar limit J-q of the 1MM free energy. From this formula and Q56JI we obtain the 
celebrated expression relating Tq and the period matrix n t j of the Riemann surface: 

' " ' dujj = Tij. (57) 



dSidSj J Bi 

Matrix integral as a Whit ham Krichever system. We now show that J-$ enjoys additional 
geometrical structures arising when differentiating w.r.t. times tk and to thus accommodating the 
whole SWWK system. 

We consider variations of the potential, i.e., variations w.r.t. (Whitham) times tk- We then obtain 
instead of (|55|) 

while still holds. Repeating the arguments of the previous paragraph, we conclude that the 
zeroes of M s -t (A) cancel from the denominator and, therefore, the variation may have pole only 
at A = oo, or r/ = 0, i.e., at the puncture. In order to find this pole, we use (|54|) . which implies 
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that dS = M m 3-t(X)y(X)dX — > (V^(A) + 0(j))dX and, therefore, the variation of dS at large A is 
completely determined by the variation of V^(A). Parameterizing V(A) = J2T=i ^fcA fc , we obtain that 

ddS , ddS ln . . 

—-=dwi, — = 2dn k 59 
dSi dt k 

with the behavior at infinity 

dn k = k (C~ k ~ 1 + 0(1)) d£, for ^ 0, k > 0. (60) 



This holds up to a linear combination of holomorphic differentials. The normalization of dVt k is fixed 
by the condition 

f)S 1 f 8S- 

^-ssV n »- 0V4 -*- - sk^- 04 - 1 (61> 

The derivatives of dS w.r.t. the times are 

o , n _ 9dS H n+k _i(X)dX 

2dtt k = — — = — — , (62) 

dt k y(A) 

and normalization conditions (|61jl together with the asymptotic expansion 

2dn k (X)\ X ->oo = kX k - 1 dX + 0(X~ 2 )dX = 

oo 

= kX k - l dX + c km X~ x - m dX. (63) 

m=l 

fixes uniquely the coefficients of the corresponding polynomials H n+k ^\ of degrees n + k — 1. 

Equations (j59j) together with normalization conditions (|61|) determine the SWWK system jjj (see 
also |4l)l |S]) whose prepotential is J-q. For this, we apply the formula similar to Q46|) allowing now 
variations of the potential V(A). We then obtain (see (|4fi )l -(|4"% |l ) 



i i 

= -res A =oo(A A ''(iS') = -v k , k = 1, . . . , m, (64) 

which is again the standard formula of the SWWK theory. 

We now consider the variation of dS w.r.t. the hard wall parameters ap. Note first that because 
Si and to are now independent variables, together with t k and an, we must demand 

SS,_ 1 I %(A) dA = ovijA md 8VW_„_ {65) 
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Performing the variation of (|24|) w.r.t. a/3, we obtain 



y(A)%y(A) = 2 , = + = = =, (66) 

Vri7=i( A - a 7) ( A -a/3)\/n7=i(A-a 7 ) (A-a^n^iCA-ay) 

so, again, using considerations on cancelation of zeros of the polynomial M(A), we have 

dy(X) P n _ x (A) 



dap (A - ap)y(\) 



dn a0 (X), (67) 



where aT2 a/3 (A) is a (noncanonical) meromorphic 1-form with the only second-order pole at A = ag and 
with all its A-cycle integrals vanishing due to condition (|65|) . To find the quadratic residue of this form 
at A = an, we need the most singular contribution when performing the variation of y(A) = M(A)y(A) 
w.r.t. an] this contribution comes from the partial derivative in an, and we have 



dn afj (X)= l -M{a p y Yla l{<1P =- dX + O (-7^=) dX; f cl9„„ A\) - V, (OS, 
2 JUU (a/3-a 7 ) (A-a/3) 3/2 VV 1 ^/ fl 

V 7#/3 



A, 



On to-dependence of the prepotential. The last property of J-q to be verified is the behavior 
w.r.t. to. From the SWWK theory, we expect that 



dt 



OCa 



dS (69) 



and = dfio where df^o is the Abelian differential of the third kind with two simple poles at two 
infinities: 

resoodfio = — reSoo_dfio = — 1- (70) 

The integral in (|69j) is however divergent and needs regularization. 

We know the derivative dTo/dto, (|49f> . which is equal to the integral § v log |A — q\dS — V(q) with 
the reference point q to be chosen on the last, nth, cut, while the expression itself does not depend 
on the actual local position of the reference point. We choose it to be «2n — the rightmost point of 
the cut. We can then invert the contour integration over the support T> to the integral along the 
contour that runs first along the upper side of the logarithmic cut from \X2n to a regularization point 
A, then over the circle C\ of large radius |A| and then back over the lower side of the logarithmic 
cut in the complex plane. In order to close the contour on the hyperelliptic Riemann surface under 
consideration, we must add the integration over the corresponding contour on the second sheet of the 
surface as shown in Fig. El we let Cl denote the completed integration contour, and it is easy to see 
that such an additional integration just double the value of the integral. 
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Figure 2: Cuts in the A-, or "eigenvalue," plane for the planar limit of 1MM. The eigenvalues are 
supposed to be located "on" the cuts. We add the logarithmic cut between two copies of the infinity 
on two sheets of the hyperelliptic Riemann surface in order to calculate the derivative w.r.t. the 
variable to and let Cl denote the corresponding integration contour. 

All the singularities appearing at the upper integration limit (i.e., at the point A) are exactly 
cancelled by the contribution coming when integrating the expression dS log(A — fi2n) along the circle 
C\; in fact, the latter can be easily done, the result is — 2ni(S(A) — (S(fi2n))+), where the function 
S(X) is the (formal) primitive of dS (which includes the logarithmic term), and we project it to the 
strictly polynomial part. Using the large-A asymptotic expansion of the differential dS, 

dS(X)\ \->oo = V'(X)dX + ^dX + 0(X~ 2 )dX, (71) 

A 

we obtain that (S(^2n))+ just cancels the term V(^2n), and we eventually find that 

= 2m I [ A dS-S(A)\ , (72) 

where Cl is the contour in Fig. [21 which by convention encircles the logarithmic cut between two 
infinities on two sheets of the Riemann surface and passes through the last, nth, cut. We have 
therefore proved that the planar limit J-q of the 1MM free energy is the SWWK prepotential, or 
semiclassical tau-function. 

We now introduce the (complete) set of canonical variables {Si, i = 1, . . . ,n — 1; to; tk, k = 
1, . . . , m}, which we uniformly denote {£/} (in what follows, Latin capitals indicate any quantity from 
this set). From (J5?i|) . (|7()|) . and (|H2*|) . we then obtain the general relation 

where Hj(X) are polynomials. 
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Asymptotic formulas (|7U|) and (|63|) imply that derivatives of all the quantities d£li w.r.t. any 
parameter are regular at infinity and may have singularities only at the branching points \ii of reduced 
Riemann surface 

Note again that these derivatives are purely geometrical and does not depend, in our setting, on 
the nature of the branching points: whether they come as hard walls or as "dynamical" branching 
points. 



2 Derivatives of the free energy and geometry 
2.1 Second derivatives of the free energy 

Bergmann bidifferential and 2-point resolvent. Calculating the one-point resolvent in the 
preceding section required the knowledge of the first derivatives of the matrix model free energy. We 
now turn to the two-point resolvents involving second derivatives. Here, instead of differentials with 
some prescribed properties of holomorphicity, the main object we need is a ^differential on a Riemann 
surface T, g : the Bergmann kernel (canonically normalized bidifferential in Fay's terminology) which 
is the double derivative of the logarithm of the Prime form E(P, Q) such that it is symmetrical in 
its arguments P,Q 6 E 9 and has the only singularity at the coinciding arguments with the behavior 
(see |H]) 

b(p,q) = ( _g(Q))2 ffC^ffCQ). ( 74 ) 

in some local coordinate £(P). As it stands, we can add to (|74|) any bilinear combination of Abelian 1- 
differentials duji] we fix the normalization claiming vanishing all the integrals over A-cycles of B(P, Q): 

I B(P,Q) = 0, for i= 1,...,<?, (75) 
JAi 

and, due to the symmetricity property, the integral may be taken over any of the variables P or Q. 

We now show that the Bergmann kernel generates the differentials aYl^- Bipole differential 1)70(1 
can be expressed through the Prime form as 

dfl = dlog — — r 76) 

E(P, oo_) 

The primitive of differential Q76JI (which we need in what follows) then obviously develops the loga- 
rithmic cut between the points of two infinities on the Riemann surface. 

Using ©, (jnU) and we obtain 

W (X,fi)d\ dfi = -Y ^^W Q ^) = ---^—V'U)dXd l i + -y ^Ay(A) 
y p; p ^ A-J+ 1 dtj VP; 2dV(X) VPy P 2^ X3 +1 dtj KJ 
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iXi, ig ^ (?7) 



2 (A- ,1)2 2^\> 

Because dQ,j(fj,) = (±j/i- J_1 + 0(l//n 2 )) dfj, for \x — ► oo±, j > 0, we conclude that the infinite sum 
3 Sj^=o jjrrd'Qj(lJ') develops the second-order pole 2 (x-^y 2 m ^ ne case wriere A an d M are on the same 
(physical) sheet; this pole is then canceled by the first term, and the second-order pole — 2 (\-fy f° r ^ 
and being on different sheets. In the latter case, this pole is doubled by the first term, and together 
with the evident normalizing conditions that follow immediately from l|77|). 

W (A, n)dfi = I W (X, n)d\ = 0, (78) 
Ai JAi 

we finally come to the formula for the 2-point resolvent, 

W (A, fi)d\dn = ^j^dXdf, = -B(P, Q*), (79) 

M^ = -l(B(P,Q)-B(P,Q*)), (80) 

where we have introduced the *-involution between two sheets of our hyperelliptic curve, so that Q* 
denotes the image of Q under this involution. The only singularity of (|79|). for a fixed point P on a 
physical sheet, is at the point Q — > P* on the unphysical sheet with fi(Q) = A(P*) = A(P). 

In what follows, we try as much as possible to express everything in terms of invariant quantities 
(the Bergmann bidifferentials and their primitives) on the reduced Riemann surface (|35[) and in terms 
of the WK differential dS. As we shall see, the 1/iV-expansion can be constructed basically in terms 
of just these characteristics. However, we somewhere need the explicit expressions for B(P,Q). 

We need a technically simpler expression for the one-differential d£(X,[i), which is the primitive 
of the Bergmann kernel B(X,n) w.r.t. the argument /i. 6 Obviously, it is a single-valued differential 
of A with zero ^4-periods on the reduced Riemann surface and is multiple-valued function of fi, which 
undergoes jumps equal to dwi(X) when the variable \i passes along the cycle Bf. 



and, correspondingly, 



d£ (^P,Q + j> J = 2TTidwj{P) + d£(P, Q), 
d£^P,Q + j^j = d£(P,Q). 



In Sec. |3 we use the explicit representation for the kernel d£(A,//) in the hyperelliptic case: 

d£ ^ a ~ (a-^(a) - g ~wr t ^w^mr (81) 



3 To be more precise, it is primitive of the antisymmetrized Bergmann kernel i(_B(A,^i) — _B(A*,/i)). 
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We take d£(X, fj,) to be antisymmetric w.r.t. the involution A —* A* between the physical and unphysical 
sheets, so it has the simple pole with the residue one at A = \i and the simple pole with the residue 
minus one at A* = // having vanishing ^4-periods in the variable A. Then d£(X,fi) = fjt Q B(X,£) — 

J^ o B(\,£), with fj,Q the reference point on which there is no dependence. 

Note that formula (|68j) expressing the derivative of y(A) w.r.t. ap is proportional to d£(X, an) and, 
moreover, can be conveniently written as the residue at the branching point ap: 

= <m ap (A) = res M=a/3 (B(X, ^)y(p)) , (82) 

whereas at a branching point fa (i = 1, . . . , 2n), we have 

2B(XM)= lim ^M. (83) 

In terms of B(P, [fii]), we have one of the Rauch variational identities |42| : 

^B(P,Q) = ~B(P,[v i ])B(Q,[tii\), i = l,...,2n. (84) 



Mixed second derivatives. Normalizing conditions for Wh(X). Another set of relations follows 
from the general properties of the Bergmann kernel, and can be also derived directly from the formulas 
of Sec.^ To this end, we apply the mixed derivatives d/dV(fi) and d/dSi to the planar limit of the free 
energy F Q . On one hand, dF /dSi = § B . dS and using that dS{\) = y{X)d\ = (V'(X) - 2W {X))dX, 

d ffv$ = — (a- m )2 an d formula (|79|). we have 

i W rf " " i ( 2B(F ' "> " ra**) = i 2B(P - ">■ (85> 

On the other hand, acting by derivatives in the opposite order, we first obtain dJ r o/dV(fj,)dfi = 
Wo(fJ,)dfi = {V'(fi) — y(fj,))dn and then d(V'(/j,) — y(fi))dfi/dSi = 2duji(n), or > i n the coordinate-free 
notation, one of the classic identities (see, e.g., [IT]): 



2m 

This means that 



— / B(P,Q) = <L)i(Q) (86) 
m JBi 

ddS(n) 



dSi 



ddS(X)_ dx 



r Bi dv^ 

where the both integrals are taken over the variable A. Now, as dS = y(X)dX is the generating function 
for the variables £ a = {{b{}, {{i a }, Ma^}, given by (O, ((23), and $Zj$ (recall that are just ith 
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order derivatives of M(A) at A = fj, a ) giving their dependence on Si and tk, one concludes that similar 
relation for the first derivatives holds for each of these variables. Indeed, multiplying (|87|) by one 
then can bring ju successively to /i Q 's, Aj's and oo to pick up pole terms with different £ a and prove 
that 

dSi = i W(x) dX (88) 

which becomes true also for any function depending on a finite number of "local" variables £ a in the 
case of arbitrary potential. For instance, it does not hold for Sj for which dSj/dSi = 8ij whereas 
8Sj/dV(X) = 0. However, one can never express Si as a function of a finite number of "local" variables 

Since, as we show in sec. El the free energy Th at any order of 1/iV-expansion depends only on a 
finite number of local variables, which are the branching points \x a and a finite number of moments 

dug 
dV{\) 



Ma \ and Wu{X) is then expressed exclusively in terms of a finite number of derivatives Jy^) and 
Wfry, we have 



(fc) 

, we have 

dT h I dTh 



dSi J Bl dv(\) 

together with the condition 



/ 7n77^ dX=( ( Wh(\)d\ Vh and i = 1, . . . ,g, (89) 
JBi oV(X) J Bl 



i ^FTrk dX = I W h (X)d\ = Mh > 1 and for i = 1, . . . ,g, (90) 
J A, dV(X) J Al 



'A dV(X) 

that follows directly from normalization condition (175(1 . 

In the same way, we can obtain the defining relation for the to-derivative 



dt Q Joo_ 8V(0 



de, (91) 



which again holds for any h and requires regularization only for h = since higher corrections in h 
are regular at infinities. 

2.2 Residue formula and WDVV equations 

The Witten-Dijkgraaf-Verlinde-Verlinde (WDVV) equations [211 HEl is the systems of algebraic equa- 
tions 

*-j -J~K = J~kJ~j 

on the third derivatives 



TtT^Tk = T K T-j X T u V I,J,K (92) 



m\\jK = 7 ^-^- = fiJK (93) 
oti at j otx 
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of some function ^({tj}). These equations can be often interpreted as associativity relations in some 
algebra for an underlying topological theory. 

That the 1MM free energy satisfy the WDVV equations was proved in (20], where it was shown that 
the multicut solution of the 1MM satisfies the WDVV equations as a function of canonical variables 
identified with the periods and residues of the generating meromorphic one-form dS [TJ. The method 
to prove it consists of two steps. The first, most difficult, step is to find the residue formula for 
third derivatives ()93|) of the 1MM free energy. Then, using an associativity, one immediately proves 
that the free energy of multi-support solution satisfies the WDVV equations if we keep the number 
of independent variables to be equal to the number of "dynamical" branching points. That is, the 
presence of the hard-wall branching points results in the reducing the total dimension of the WDVV 
system. 



Residue formula Because all the quantities d£li ()73l) depend entirely on reduced hyperelliptic 
Riemann surface 1)35 Jl. their derivatives w.r.t. any parameter must be expressed through the derivatives 
w.r.t. the positions of the branching points We skip most of details referring the reader to [71 I2U|. 

The formula for the third derivative d s J r o/(dtidtjdiJ a ) = Fjja easily follows from (|8l|) when tj 
and tj are times of the potential, and one should use the Riemann bilinear identities in the case of 
other canonical variables. 



We eventually have 



^^(W;) ^^, (<*) 



Completing the calculation of the third derivative needs just inverting the dependence on the 
branching points therefore finding d[i a jdtK- (Note that, obviously, dap/dtx = 0.) Differentiating 
expressions (|24|). ()25|) w.r.t. tpc, we obtain for (|59|) 

ddS(X) _ H K (X)dX 
dt K ~ y(A) 

= ^(A) ± jf^pdX + »(A)dA. (95) 

2 a^l ( A ~ Vol) Ot K dt K 

The derivative of the polynomial M(X) is obviously polynomial and regular at A = fj, a . Multiplying 
(|95|) by y/X — [x a and setting A = fj, a , we immediately obtain 

djM = H K (fi a ) 

m K M(/i a )nw(/i«-M7)' 

with only "dynamical" branching points in the product in the denominator. Combining this with (|94j) . 
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we come to the desired residue formula for the third derivative w.r.t. the canonical variables ti: 



d 3 Fo _ Hl{Ha)Hj{n a )H K {ii c 



dt I dt J dt K ^ M(p a ) n S 7=i (/*<* ~ ^ 7 ) 2 IljUiCMa - a p) 

•J-^ dQjdUjd^lK dQidUjdQx dQidQjdQ,K 

= X>*. d Xdy = dAdy =reSdA =° dAd, • (97) 

We are able to pass to summation over all branching points in the last line because the residues at 
A = ap just vanishes and in the very last expression we assume zeros of the differential d\ on the 
reduced Riemann surface (|35|) . 



Proof of WDVV equations Given residue formula (|97[). the proof of WDVV equations (|92|) 
can be done by checking associativity of the algebra of differentials dVti with multiplication modulo 
x/ria=i(^ ~~ n a )d\. This algebra is reduced to the algebra of polynomials Hj(X) with multiplication 
modulo []q=i(A — MoDj which is correctly defined and associative. The basis of the algebra of -ffr(A) 
obviously has dimension s and is given, e.g., by monomials of the degrees 0, 1, . . . , s — 2, s — 1. 

Another proof pertains to solving the system of linear equations [431 144j . Defining 



we reduce (|97j) to the form 



F?)k = E (99) 

a=l 



We also demand the determinant of cf>J to be nonzero: 

det||0?|| #0. (100) 

la 

This nondegeneracy holds automatically in the case where we have no hard walls and choose the 
canonical variables to be Si, io 5 an d first n times (k = 1, . . . , n) of the potential [U- 7 In the case 

7 From the definition 1981 . this condition stems, up to obviously nonvanishing Vandermonde determinant factors, to 
the nondegeneracy of the matrix a 11451 1. Indeed, if deter vanishes, then there exists a polynomial P(A) of degree less or 
equal n — 2 such that 

• P(A) 



A»2i-1 av ' 



d\ = for i = 1, 



This necessarily implies that P(A) has at least one zero at each of the intervals (fl2i-i, 1^2%)', otherwise the combination 
under the integral sign is sign definite and the integral cannot vanish. The polynomial P(A) must then have at least 
n — 1 zero and, having the degree not exceeding n — 2, must therefore vanish. 



23 



with hard walls, however, no such consideration works, and this determinant may vanish on subdo- 
mains of parameters of the problem, so we must consider condition (|lUUj) as an additional restriction, 
besides the matching condition claiming that the number of canonical variables must coincide with the 
number of "dynamical" branching points of the solution. Recall that we need at least three canonical 
variables in order to have a nontrivial WDVV system. 

We now let J-f 1 ^ denote the matrix with the entry JK equal to J'fjx- The WDVV equations Q92j) 
follow now immediately from (|99|) and nondegeneracy condition UlUUj) , 

3 Higher genus contributions 

The solution Wi(A) to the loop equations in the multicut case was first found by Akemann |35j and 
the universality of critical behavior of the corresponding correlation functions was demonstrated in 
Akemann also managed to integrate the 1-point resolvent to obtain the free energy T\ in the two-cut 
case. The genus-one partition function in the generic multi-cut case was proposed in |45| I4fi| . where 
it was observed that the Akemann formula coincides with the correlator of twist fields, computed by 
Al.Zamolodchikov j^Jj. The function F\ in the general multicut case was found by integrating Wi(A) 
in @H], HQ. 

3.1 The iterative solution of the loop equation 

Inverting the integral operator. We now determine higher genus contributions. We do this 
iteratively by inverting the genus expanded loop equation (|15|). Our strategy will be to construct an 
integral operator d£ inverse to the integral operator K — 2Wo(A). 

Acting with this operator on the both sides of loop equation (|15|) we recursively obtain Wh(X) for 
all genera like all the multi-point resolvents of the same genus can be obtained from Wh(X) merely by 
applying the loop insertion operator gy^ ■ Another ingredient of the construction will be the action 
of the loop insertion operator within our technique. 

However, there is a subtlety: the operator K — 2 Wo (A) has zero modes and is not invertible. 
Therefore, solution to the loop equation is determined up to an arbitrary combination of these zero 
modes. Moreover, in the case of hard walls, the kernel of d£ does not coincide with the one of 
K — 2 Wo (A) and must be discussed separately. 

From representation (|18|) and l)23j) . we have that the action of K — 2Wq on W/ l (A) is 



(K-2W (X))W h (X) 



[V'(X)W h (X)]_-(V(X) - y(X))W h (X) 
[V'(X)W h (X)U- [(V'(X)-y(X))W h (X)] 
[y(X)W h (X)]_. 



(101) 
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Because both Wh{X) and y(A) are odd functions under the involution A —* A* interchanging the 
physical and unphysical sheets, we conclude that the combination y(X)Wh(X), being an even function 
under this involution, is a meromorphic function on the physical sheet, that is, it is a rational function 
having poles only at the branching points and at infinity. 

It is now easy to see that the integral operator 

with the kernel d£(X,/j>) from (|81jl is an inverse for the operator K — 2Wq{\). Indeed, representing 
f(p) in form (jlOlj) . we have 



A>c D 27ri dX y(ji) T\>Ct, 2tti dX y(fj) 

dfj. d£(X,fi) ( w /^ Pm-2(j*) 



\ >Cv 2ni dX V y(fi) 

Jc 00+ 27rz dX J\>Ct, 2m dX y(fj,) 

where P m _2(/i) is a polynomial. Here, we indicated by A > Cx> that the point A lies outside the 
integration contour in the physical sheet and, when pulling the contour of integration to 00+ , only the 
residue at /x = A contributes in the first integral and, by the normalization conditions, this residue is 
exactly Wh(X). The remaining integral at infinity vanishes because the integrand is regular as \i — > 00. 
Next, by virtue of explicit formula (|81|) . we see that because the last term is a rational function on 
the physical sheet that is regular at the branching points, integrating it over Cv obviously gives zero, 
so the whole action of d£ on (K — 2Wo)W/ l (A) just reconstructs Wh(X). 

Note also that under the action on [y([i)Wh(fJ-)]-, due to the involution properties, we can replace 
the integration over Ct> just by evaluating residues at fj, = m in formula ()1(J2|) having 

d£(f)(X)dX=i ^^(A^-^./Ox), (104) 

where we let denote the union of contours encircling all the branching points. 
Because the operator d£ obeys the property 



$ d£(f)(X)dX = (105) 



it also enjoys integrability conditions (|90j) . Therefore, one has to solve the loop equations inverting 
K — 2Wq exactly with d£. 
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Zero modes. The modified loop equation. Note that now, in contrast to the case without hard 
walls, the zero-mode content of the operators K — 2Wq and d£ is different. 

It is clear from representation (jlUlj) that the zero modes of the operator K — 2Wq when acting on 
the functions Wh(X) with asymptotic behavior (1161) are spanned by P s -t 9 (X)/y(X) and constitute the 

2 

linear space of dimension — 1. In addition, all the simple-pole functions proportional to x-ap are 

annihilated by the operator d£. Altogether, it gives (in the case s > t) the total space of zero modes 
of dimension — 1 + i = re — 1 = g, and the total number of zero modes under the action of the 
combination of the operators d£ o(K — 2Wq) is therefore equal to the number g of Abelian differentials, 
whereas the freedom to add these differentials is completely fixed by normalization condition (JUUJ). 

A crucial observation is that when acting by the operator d£ on the both sides of loop equation 
1)15 Jl . the terms with the partial derivatives of Fh w.r.t. ap vanish, the action of K — 2Wq on Wh(X) is 
inverted, and the remained expression just becomes 

/h-l o \ 

W h (X) = d£ \J2 W h *(\)W h - h ,(A) + w(XjW h -i(\)j , (106) 

with Wh{\) now automatically satisfying normalization conditions ()90(l . 

It follows from the properties of the operator d£ that, for all s > at h > and for s > 3 at h = 0, 
the s-point resolvent W/ l (Ai, . . . , A s ) possesses the gradation w.r.t. involution between variables of the 
physical and unphysical sheets: 

W h {\i, ...,Xj,...X s ) = -W h (Xi, ... ,X*,. . . X s ). (107) 

In particular, Wh(X, X) as well as any product Wh(X)Wh'(X) in the r.h.s. of (|106[) are then rational 
functions of A for any h, h! > 0. Moreover, Wo(A, A) = B(X, A*) is again a rational function of A having 
poles of the second order at the point A = /ij of merging of the two sheets. 

However, in order to have loop equation (|15|l. not (|106|) . we must demand additional conditions to 
be satisfied by the matrix-model free energy. In contrast to conditions (|89j) satisfied by any function of 
local variables, these conditions impose restrictions; indeed, it is required for the parts of loop equation 
(j!5l) proportional to x-ap *° ma tch. That is, the constructed loop means Wh(X) must satisfy 



da 



£ 2^ (~f + E w h ,(0w h - h >(0 + a^Wh-i(o) , (ios) 



a /3 \ h'=l 

because the residue at ^ = ap of the positive part [y{€)Wh{Q] + is obviously zero. 
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In two next paragraphs, we solve recurrently modified loop equation ()1(J6|) . We then integrate it 
to produce the free energy term T\ t in the next section and, after obtaining the answer, must verify 
conditions (|1(J8|) . Note that the two normalization conditions, (|89[) and (|ll)8j) . leave no room for adding 
terms independent on V(£) and dependent either on Si or on an alone. And, as we show in the last 
section, in the special case of T\, condition (|1U8|) fixes the function of ap to be added to the free 
energy. 

Diagrammatic technique. The above considerations provide a basis for the diagrammatic rep- 
resentation for resolvents, which almost literally repeats the diagrammatic technique for the 1MM 
without hard walls }22j . Let us represent the one-form d£(p,q) as the vector directed from p to q, 
the three-point vertex as the dot in which we assume the integration over q, • = § ^ 2y(q) > an< ^ ^ e 
Bergmann bidifferential B{p,q) as a nonarrowed edge connecting points p and q. The graphic repre- 
sentation for a solution of (|1U6|) then looks as follows. Representing the multiresolvent Wh'ipi, ■ ■ ■ ,p s ) 
as the block with s external legs and with the index h! , we obtain [22] 




(109) 



which provides the basis for the diagrammatic representation for Wh(pi, ■ ■ ■ ,p s )- It can be formulated 
as a set of the following diagrammatic rules |22j . 

The multiresolvent Whipi, . . . ,p s ) is presented as a finite sum of all possible connected graphs with 
h loops and s external legs such that 

• only three- valent internal vertices are allowed (the total number of edges is then 2s + 3h — 3, and 
we assume s > 1 for h > 1 and s > 3 for h = 0); all the vertices are indexed by the corresponding 
coordinate variables q on the Riemann surface y(A). 

• We have propagators (edges) which either connect internal vertices (can start and terminate at 
the same vertex) or come as an external leg to an internal vertex. 

• These propagators are of two sorts: there are arrowed propagators and nonarrowed propagators. 
As above, we set the 1-forms d£(p, q) into the correspondence to the arrowed propagators and the 
variables p and q correspond to the respective vertices at which this edge starts and terminates. 
We have the Bergmann bidifferential B(p, q) corresponding to each nonarrowed edge that starts 
and terminates at the vertices labeled by the respective variables p and q and if it is the same 
vertex, then B(p,p*) corresponds to this propagator. 
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• The arrangement of propagators is subject to restrictions. We have exactly 2k + s — 2 arrowed 
edges; these edges must constitute the maximum rooted tree subgraph with all arrows directed 
from the root. Different choices of the tree subgraph for the same graph are considered different 
cases of the construction, and we must consider all possible cases. The root of the tree subgraph 
is always a selected (one and for all) external leg, say, p\ (the choice is arbitrary due to the 
symmetry of Wk[pi, ■ ■ ■ ,Ps)) and the maximality of the tree subgraph means that it contains no 
loops (including tadpoles; the arrowed propagators cannot therefore start and terminate at the 
same vertex), is connected, and for every internal vertex there is exactly one arrowed propagator 
that terminates at this vertex. We call this propagator the incoming propagator for this vertex. 
There can be either one, or two, or no outgoing arrowed propagators at an internal vertex of the 
graph. 

• Each tree subgraph establishes the relation of partial ordering on the set of internal vertices 
of the graph; among two vertices one precedes to the other if there exist a naturally oriented 
path in the rooted arrowed tree subgraph that starts at the first vertex and terminates at the 
second vertex. At each internal vertex we place the integration § c { q ) ^2y(gJ over ^ ne variable 
indexing this vertex; the integration goes along the contour that is about the domain T> and 
the arrangement of the integration contours at different vertices is prescribed by the arrowed 
subtree: the closer is a vertex to the root, the more outer is the integration contour. 

• All internal nonarrowed propagators are allowed to connect only vertices that are in the partial 
ordering relation to one another; it can be the same vertex. All external propagators except 
the one that is the root of the tree subgraph are nonarrowed, and we have the Bergmann 
bidifferentials B(q,pk) (k = 2, . . . , s) corresponding to these edges. 

• All the integration contours over the A-cycles in formula (|81jl for d£(X,fj,) are assumed to be 
outside all the integration contours over the variables indexing the internal vertices whereas all 
the external points pi, i = 1, . . . , s are assumed to be outside all the integration contours C$ 
(as well as outside all the contours of integration over the ^-cycles). 

We now demonstrate the consistency of this diagrammatic technique by calculating the action of 
loop insertion operator d/dV(r) on its elements following |23j . 

We first calculate the action of d/dV(r) on B(P,Q) thus producing an analogue of formula (|97|). 
Using (|84j) . we represent this action through the action of partial derivatives in fi a (a = 1, . . . , s; only 
the dynamical branching points participate on this stage) subsequently calculating the latter from 
relation (jHOJ- Let 



y(x)dx\ x ^ a = y(\p a ])-\/x - \x a dx + 0(y/x - n a ) dx 



Then, since 



dy(p)dp 



dV{r) 



2"^ 1J y/p-=J^dV{r) 
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we have 

d[x a _ 2B([fj, a ],r) 
dV(r) y([» a ]), { > 

and, therefore 

9V(r) ^ V(M) 

This is apparently just formula (|97jl upon setting all variables /, J, K to be the times of the potential 
with subsequent summation with p~*~ 1 , q~^~^, and r~ fc_1 for p, q, and r being the variables on the 
physical sheet. 

However, now we represent expression in the different way. Instead of differentiating y(£) in 

the denominator as we did in formula (|97|1. we integrate one of the Bergmann kernels, which gives the 
1-differential d£(Q,£). Note that the local variable in the vicinity of fi a is £(x) = sfx — (j, a , this gives 
the additional factor 1/2, so we have 

-B(P,Q) = 2^res MQ - 



9V(r) v ^ ^ 2y(i{x))dt{x) 



and we can again add residues at the hard wall points an because they are just vanish here, so our 
final expression is 

8 2n 



dV(r) v fri ^ 2y(£(x))d£(x) 

The price for preserving y(£(x))d£(x) = dS in the denominator is that we have lost the explicit 
permutational symmetry for the three terms in the numerator, which was apparent in (|97|), We 
however pay this price for the possibility to develop the comprehensive diagrammatic technique. From 
this relation, it obviously follows that 

J-m Q , P) - f res, ■*<«»), P »t W £f B(r - «'» , (113) 

and the last quantity to evaluate is 

d 1 B(p,r) 
dV(r)2y(p) = ~2y 2 (p)' 



(114) 



Note that the point P in (|113|) is outside the integration contour. Multiplying the both sides 
of (|113|) by l/(2y(P)), using (|114|) . and pushing the integration contour through the point P, we 
observe that the contribution of the simple pole of d£(£(x), P) at the point £(x) = P cancels exactly 
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the variation of 1/ '(2y(P))! We therefore attain the prescribed contour ordering and can graphically 
present the action of d/dV(r) as 

r r r 

^4) = Q^L^, m%) <3 p = Q -J— p = Q—^P- ( 115 ) 

In the second case, we have the freedom to choose on which of edges to set the arrow. Recall, however, 
that as we have had a nonarrowed propagator stretched between the points P and Q, these points 
were to be ordered before, as prescribed by the diagrammatic technique. That is, if "P preceded Q" 
we must choose the first variant and if U Q preceded P" we must choose the second variant of arrows 
arrangement. 



3.2 Inverting the loop insertion operator. Free energy 

The iif- operator. We now introduce the operator that is inverse to loop insertion operator lj7jl. 
Let 8 

1 1 r°°+ n_1 r 

H- = -Te Soo+ V(x) ■ --ve Soo _V(x) ■ -t / • - V S t <b • . (116) 

I 2 Joo_ J Bi 

The arrangement of the integration contours is as in Fig. We calculate the action of H on the 
Bergmann bidifferential B(x,q) using again the Riemann bilinear identities. We first note that as 
B(x,q) = d x d£(q,x), we can evaluate residues at infinities by parts. Then, since d£(q,x) is regular at 
infinities, we substitute 2y(x) + 2to/x for V'(x) as x — ► oo + and —2y{x) + 2to/x for V'(x) as x — > oo_ 
thus obtaining 

-reSoo + ^y(x) + ^ (g, x)dx + reSoo_ y— y{x) + d£(q, x)dx 

n—l „ 

-t d£(q,x) -J^Sif B(q,x). (117) 

x= °°- f^l jB i 

Whereas the cancelation of terms containing to is obvious, it remains only to take the combination 
of residues at infinities involving y(x). For this, we cut the surface along A- and P-cycles taking 
into account the residue at x = q. The boundary integrals on two sides of the cut at Bi then differ 
by d£(q,x) — d£(q,x + § A ) = 0, while the integrals on the two sides of the cut at Ai differ by 
d£(q, x) — d£(q, x + <f B ) = § B B(q, x), and the boundary term therefore becomes 

n-1 



f y(x)dx j> B(q,£), 



This definition works well when acting on 1-forms regular at infinities. Otherwise (say, in the case of Wo(p)), the 
integral in the third term must be regularized, e.g., by replacing it by the contour integral around the logarithmic cut 
stretched between two infinities as was done in Sec. 11.51 
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which exactly cancels the last term in pi7|) . Only the contribution from the pole at x = q then 
survives, and this contribution is just —y(q). We have therefore proved that 



H-B{;q) = -y{q)dq. (118) 

Free energy. We now act by H on W&(-) subsequently evaluating the action of loop insertion 
operator (JJJ) on the result. Note first that the only result of action of d/dV{p) on the operator H 
itself are derivatives dV(x)/dV(p) = — l/(p — x) (and recall that by definition \p\ > \x\, i.e., instead 
of evaluating residues at infinities one should take residues at x = p) , which gives 

(H ■ W h (.)) = W h (p) + H ■ W h (-,p). (119) 

For the second term, due to the symmetry of Wh(p, q), we may choose now the point p as the root of 
all the tree subgraphs. Then, the operator H always acts on £>(•,£) where £; are integration variables 
of internal vertices. However, if this vertex is an innermost (i.e., it has no outgoing arrowed edges), 
then the 1-form y '{Cjd^ arising under the action of H (|118|) cancels the corresponding 1-form in the 
denominator of the integrand, which becomes therefore regular at the branching point giving zero 
contribution. If this vertex has an outgoing arrowed edge, say d£(£, p) (it can be at most one outgoing 
arrowed edge as the third edge must be external), then, again, we can push the integration contour 
for £ through the one for p; the only contribution comes only from the pole at £ = p. The value of 
the residue is however doubled, and we come to the following graphical representation for the action 
of the operator H: 



H- H- 

Q-J-^p = - Q — -p; Q-J— p = o. 



(120) 



For H q ■ Wk(q,p) = H q ■ gy/ q \ Wk{p), we obtain that for each arrowed edge on which the action of 
d/dV(r) produces the new vertex, the inverse action of H q - just give the factor —1 and for each 
nonarrowed edge on which the action of d/dV(r) produces the new vertex, the inverse action of H q - 
just gives zero. As the total number of arrowed edges is exactly 2h — 1 for every graph contributing 
to the sum of diagrams with one external edge and containing h loops, we obtain that 

H q -W h (q,p) = -(2h-l)W h (p). 

Therefore, combining with ()119j) . we just obtain 

(H q -W h (q)) = -(2-2h)——F h , (121) 



dV{p) K q nvi " v 'dV{p)' 
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and, since all the dependence on filling fractions and to is fixed by the claim that the answer depends 
only on ^ and derivatives of M(p) of finite orders at the branching points, we conclude that 



* = 2h=2 H ' Wh - (122) 

We have therefore expressed the free energy through the one- loop resolvent. Using this final answer, 
we can calculate all Th except the contribution at h = 1 (the torus approximation). In the case 
without hard walls, the latter was calculated by a direct integration in |48j (see the refined version 
in [25)) an d we generalize this calculation to the hard wall case in Sec. HI All other orders can be 
calculated consistently. For this, we need to introduce one new vertex © in which we place the new 
integration term: 



i Mi J 



(123) 



This is because, although the integral term Jn o y{s)ds is nonlocal, its constant part Jq q , having the 
improper involution symmetry, drops out of the residue in the 1MM case, and we can integrate only 
in the neighborhood of each branching point \ii separately. That is, we need only the local expansion 
of this integral as £ — > /^. Then, say, the genus two contribution is provided by the sum of three 
diagrams 



2 .^ 2 = 2 o*o +2 (V) +CHH3- 



(124) 



Note also that J-q written in form (|52|) obviously can be presented as 

Fo = --H xee -W (125) 

due to (|48|) and l)49j) and by virtue of relations (|89|1 and (j91j) . Recall that in the case of J-q, the integral 
between oo_ and oo+ must be regularized as in Sec. 11.51 

Postponing the consideration of the case T\ to the next section, let us show here that H ■ W\ is 
constant. Indeed, 

^ In V(s)ds 

~i viz) 

and the first term has simple zero with residue 2/3 or 2 (depending on which branching point — hard 
wall or dynamical — we have at /ij) to be compensated by the double pole of !?(£,£*) ~ 4^-^ p ' an< ^ 
the total answer is then just the constant equal to (2s/3 + 2t)/4. 
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Scaling relation. We now explain the origin of relation 1)121(1 , Indeed, recall that for any functional 
T of a finite number of "local" variables, which are in our case the branching points pi and the 



moments , we have relations ()89|) and ()91|) that state that 



1 



00 -1 

OO — 



d£ = dJ-h/dtQ. Relation (|121l) is then equivalent to the formula 



x 8 8 n ~~ 1 8 8 



'dh 



dTh/dSi and that 



(126) 



which is just a ("kinematic") property of every integral of form Indeed, let us fix the total number 
of eigenvalues to be N and the partial numbers of eigenvalues on intervals of distribution to be iVj. 
Then, to = hN, Si = fiNi, and we have the exponential term \ Y^k^k^ X k . If we just keep N, 2Vj, 
and all the combinations t^jfi fixed and scale only the formal expansion parameter h, we immediately 
come to (HSU). 

Note that in the case of 1MM without hard walls, the operator in (|126j) has also the sense of the 
scaling transformations under which all vertices are multiplied by the same scaling factor p and all 
the propagators are multiplied by p . The action of derivatives in to an d Si results in multiplying 
all index loops (faces of the fat graph) by p. Therefore, for any graph, the total factor is 



# vertices - # edges 



# faces ^2- 



2h 



and it is exactly canceled by the scaling of the formal expansion parameter fi — > ph. 

The above relation, being purely kinematic, pertains to any matrix integral irrespectively to the 
domain of eigenvalues, etc.; to the best of our knowledge, this is the first time when it found a technical 
application. 



On checking the consistency conditions for h > 2. The consistency condition for h = are 
satisfied by construction. We verify these conditions in the case h = 1 in the next section and consider 
here only the case h > 2 in which we have now the explicit representation for the free energy Th- 

A simple, but rather naive, consideration pertains to that the whole expression for J~h with h > 2 
can be treated as a functional of y(p). Then, 

/ ^b/(A)Wfc(A)] 



dap 



f dX_ dT h f dX f dp dy{p) SJ^h 

J Cafi 27ti y[ dV(\) ~ Jc ap 2^I x>Ct> 27ri y[ dV(X) 5y(p) 

r dp dy{p) bT h = f d\_ f dp, A)y(A) — (128) 

7c D 2ni dap 5y(p) Jc ap 2iri 7a<c c 2-iri ' 5y(p) ' 
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and these expressions differ only by the order of integration. 9 
Let us introduce the operator 








(129) 



dV{\) ' 



When acting on functionals of y(/x), it can be presented as 




a<Ce, 27ri 





-^B( f i,X)y(X) 



(130) 



and we can schematically present the action of Dp as 





(131) 



where we let the white dot represent the function y(fj,) and the fat contour designate the integration 
over the circle about the branching point ap; apparently, integrations about other branching points 
commute with this integration. We can now use representation (|106ft . which implies that, actually, 
when acting on Th, we can put all variations to the outer line to obtain, instead of (|131|) . 



where the asterisk now designates the whole combination Yjy=\ W^(X)Wh-h'{X) + Wh-i{X, A) in the 
r.h.s. of (|l()(i|h However, for these considerations to work we must be able to collapse integration 
contours to the branching points to produce representation ()1U9|) for the one-loop resolvent. So the 
above considerations are only plausible reasonings, and in order to prove the consistency conditions 
rigorously we must use other technique. At present, the only systematic way to prove these conditions 
is to use the induction in the framework of the above diagrammatic technique. This proof is however 
too cumbersome to be presented here and will be published elsewhere. 

9 We also can substitute just y{\)Wh{^) for [y(\)Wh(ty]- when evaluating the residue at A = ap. 





(132) 
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4 Calculation in genus one 



This is the only place where we integrate the expressions for Wh{\) explicitly in the spirit of [31] and 

ESI- 

4.1 Finding the free energy 

Choosing a special basis. We now express the action of operator (|l(J4j) on the monomials (/i — /ij) 
in terms of the derivatives w.r.t. d/dV(fj,). 

Using the above conditions and formula (|9(J|) . we can invert the operator K — 2Wq{\) when acting 
on basis monomials (A — fj, a ) . We define the basis vectors x\ W to be 

X ( t ) (X)d\=<f ^-^r^-r^-r^rsUX-^), i = l,...,2n. (133) 

Since the r.h.s. of (jl()6[) is this operator acting on a collection of monomials, the general form of the 
one-point resolvent is 

3/i-l 2n 

W h {\) = £ HA ( Sxf°(A), h > 1, (134) 

k=l 1=1 

where A^l are certain functions of fi a , ap, and the moments M- k \ As the order of the highest 
singularity term 1/((A — /i a ) 3/l ~ 1 y(A)) in Wh(X) is insensitive to a multi-cut structure, Wh(X) will 
depend on at most 2n(3/i — 2) moments, just like the one-cut solution ^33] • Here we need only few 
first basis functions. 

To obtain Xa\ let us start with (|%U|) . It implies 



1 d[l a 



A - /j, a dV(fi) 



= 2(S(A, M ) - B(\*,ri)-±- +0(1), (135) 



whence, as / 



2VA- 










A— » 



1 d£(iJ,,X) 1 



. (136) 



>c„ a 27ri A-// Q y(A) VA-/XC 
It then technically useful to introduce the function 

y<x(v) = v{y){^ - Va)~ 1/2 , (137) 
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which, together with all its derivatives in fj,, is nonsingular as fi —* fi a . We then apply integration by 
parts to obtain 



dfi d£(A,(j,) 



1 



dfi 1 d /d£(A,/i) 
C { } 2m y(/i) (fi-fi a ) 2 J Ctia 2iri fi- fi a dfi\ y{fi) 



dfi 1 f(B(X,n)-B(X*,f Jl )) d£{X,fi) 



C MQ 2vri fi- fi a \ 



+ 



i 



y Q (/i) 2(/j-// a ) 



(138) 



Taking the term with fi — fi a in the square brackets in the r.h.s. to the l.h.s., taking into account that 
the term with the Bergmann kernel in the r.h.s. is just —dy(fi)/dV(X), and using (|136j) . we obtain 



dfi d£(X,fi) 



1 



dpi 1 



2 -^W 2vri y(fx) (fi - fi a ) 2 J Clla 2iri fi - fi a \ dV(X) y(fi) y(fi) 



dy(fi) 1 + d£(X,fi) 



d_ 

dfi 



log y a (v) 



dy 1 



c MQ 2vrz /x - fi a \ dV(X) 







log y a (fi) + 



1 



dfi. 







C MQ 2vri V M - Ma #V(A) 
d/i 1 



(A) ^/c^ 2vri /i - fi c 
d 

(log yaWa)) 



log 2/«(/i) 



logy«(M) 



2(/i - /J 

dfl a 



dfl a d 



dV{X) (/i - ^) 2 



logy a (/i) 



log J/ q (m) 



ay(A) 

The last nontrivial action to consider is 

dfi d£(X,fi) 



On the other hand, 




dV{X) 



M{ap) 



dy y(jj) 1 



t B{\,[ap\) 







9V(X) J Ca/3 2m fi-ap y{fi) 

dfi f-B(X,fi) 1 M{fi) 



+ 



E 



c a/3 n- a p y(p) fi-ap^ 2(p - fi a ) dV(X) J 



B(X, [a /3 ])M(a /3 ) 



y(M) 



+%)E 



<9/j c 



i 2(ap-n a ) OV(X) 



Collecting all the results, we have 



1 

X- Ho 



dV(X) 



i & I? • • • > 



(139) 



(140) 



(141) 
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2 5 / 


(A 




3dV(\) \ 








A - 


1 


13 = 1,... 








-2 ° ( 


(A 


-^) 2 ) " 


dV(X) { 



11/3=1 v 7 /^ - ap 



4 2) ( A ) = ^ ( a -„ y» ) = -^K^' II , /?= I,---,*. (142) 

a=l 



Integrating Wi(X). Now we evaluate the action of dS on Wo (A, A) in the special basis above to 
obtain the subleading term J-\ for the free energy. We first obtain the explicit formula for Wo(A, A) by 
differentiating expression (|8Tj) w.r.t. fi, subtracting the term 1/(A — fi) 2 , and taking the limit [i — ► A: 



dv(\) uw 4y(A) 4^ tv y IftA-^-J^ s (e-/ij)i/(0 



1 2n 1 1 1 / 1 1 



1 ^ 1 1 ^ 1_ 

ifi 2^ M _ „.\2 Q 2^ ,,. _ 



16^(A-/Xi) 2 8 ^ fa - nj \\ - m \ — Hj 

1 n— ^ { 2n \ f I 



4S^ A Hgx^i^ra(i))' (143> 



Turning to the last term containing hyperelliptic integrals, we now prove that 



2^-logdeta, j = l,...,2n, (144) 



where the matrix a has entries 



(Ti,k=f ^ttt^j ?,fc = l,...,n- 1. (145) 

y(£) 

For the canonical polynomials -Hfc(£) = StTi -^,fc^ -1 > /c = 1, . . . , n — 1, related to the canonically 

W 

n-1 

X! a i,i H i,k = $i,k for i, A; = 1, . . . , n - 1. (146) 



normalized differentials (|56j) . i.e., such that ^ffi ^£ = &k,i-> we obviously have 
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Therefore, for all k > such that I — k — 1 > 0, 



n-l 
i=l 



and 



n-l 



n— 1 „ *r^n— 1 tt v^i— 1 ti— fe — l,,fc— 1 



(147) 



so we obtain the first part of 1)144(1 . To obtain the second part, note that by virtue of 1)146)1 it is just 

n-l 







so we come to the second equality in (|144[) . 

To find W\(\) it remains to act by d£ on 1)143)1 using 1)142)1 . This eventually gives 



Wi(A) 







s t 

.-io g nM«--iognM«--iog n 



a=l 



(3=1 



dV{\) 

~Q iog n n ~ ~ o io g det ° 



l<a<7<s 



(148) 



Q = l 0=1 

This expression is easy to integrate. Up to a function T{\a^\) depending only on a/3, we have 



t x = -i log f n M0i„) • n ^ 3 m • a(^) 4 • n(^ - ^ 

y a =l 0=1 a ,f3 



( det it 

i,j=l,...,n-l 



\12 



(149) 



where A(/x) = IIi<a<7<s(^a — M7) i s the Vandermonde determinant. This is our final answer for the 
genus-one partition function. A natural candidate for a function ^-({ap}) in this case might be a power 
of the Vandermonde determinant composed from an. But, as we demonstrate in the next subsection, 
explicitly solving the consistency conditions, we obtain ^({ag}) = and so result 1)149)) is exact. 

4.2 Verifying the consistency conditions 

Using operator Dp 1)129)1 . we can present consistency conditions 1)108)) in the case of T\ in the form 



(150) 



38 



with T\ from ()149() and Wo(£,£) from ()143|) . To obtain the answer we need only to define the action 
of Dp on the moments and on the branching point parameters. First, from (|82jl. taking the leading 
term of B(\,/j.) as A — > fj, a , we immediately obtain 



Dp/J, a = 0, a=l,...,s, 



Dm 



l,...,t. 



(151) 



(k) (k) 

It is also easy to obtain the action of Dp on moments M a and M) for 7^/3. Representing 



(AO 



C M , 2vri (^-/ii) fe y(0' 



in the case fa 7^ 0/3 we obtain that Dp gives zero when acting on all fj, a in the denominator of the 
integrand, the contours of integration about ap and fa do not intersect, so the action of Dp on y(£) 
also vanishes, and the only contribution comes from the explicit variation of just the parameter ap, 
which eventually gives 



D R M, 



(fe) 



1 1 



d*- 1 ( M(g) \ 



2 {k- iy.de- 1 \i-ap) 



fa / ap. 



(152) 



(k) 

The last quantity to calculate is the variation of M\ . Here, as above, all variations in /i a disappear, 
but contribution when varying is nonzero because of contour ordering. Permuting contours of 
integration as in (|131[) . we obtain extra derivative in £, so, denoting by 5 a „ the explicit variation w.r.t. 
the parameter ap, we obtain 



1 



y(0 



d 1 



(£ " ap) k m Jc ag 2m (£ - apf-^ d£ ^(C) ' 



where, as in (|137|) . we let yp(£.) — — a p denote the regular part of the expression y(£) in the 

local coordinates near ap. We then eventually obtain 



DpM^ 



Qk- 



( 



(k - 1)! d^' 1 



1 * M(g) 1 * M(0 



7 =i £ - a 7 



/3=1, 



(153) 



Substituting these formulas in Q15UJI and evaluating the both sides of the equation using explicit 
formulas Q149J) and (|143l) . we find the exact coincidence, which means that no additional factors 
depending exclusively on the hard wall parameters ap are required to satisfy the consistency conditions. 
Answer ()149J1 is therefore exact, and we show below that it also agrees with the ideology of Bonnet, 
David, and Eynard jH]. 
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4.3 Genus-one free energy, determinant representation, and singularities 



We now discuss the result obtained for the genus-one free energy. Comparing it with the "old" matrix- 
model two-cut solution |35| we find the coincidence up to the modular transformation that permutes 
A- and l?-periods. This is exactly because we set all the ^-periods of dS constant under 
the action of gy^ , see (|61[). On the contrary, if we impose the condition of zero -B-periods of dS, 
which corresponds to the equal "chemical potentials" IT in different wells of the potential |49| . or 
to additional minimization of free energy IJ43JI w.r.t. occupation numbers (|45)) . 10 we come to the 
modularly transformed basis of cycles. In this case, the matrix of the ^4-periods of must 

be replaced by the matrix of the corresponding .B-periods ^7j- Certainly, (|149|) reproduces the 
answer of |451 146| n for the generic multi-cut solution. 

That the only result of interchanging A- and -B-periods is the interchanging of the corresponding 
periods in a implies that e^ 1 is a density, not a function, on the moduli space of curves. Indeed, the 
corresponding determinants are related as detcr^ = det r det a^ A \ where Tij = § B . dvjj is the period 
matrix, which is itself related to J-q (see (|57|0 . In order to compensate this determinant e^ 1 must 

1 /2 

be transformed under interchanging A- and i?-cycles with the additional factor (detrjj) ' , which is 
character for a section of determinant bundle of the operator 3 over the moduli space. 

It was proposed in [12] that, in order to match the proper behavior under modular transformations, 
the operator 8 must act on the twisted bosons on the hyperelliptic curves; the factor e^ 1 is then equal 
to its determinant. Besides, one also needs to add some corrections due to the star operators [HjI I5U] 
that do not contain det a factors and cannot be restored by the modular invariance alone. These 
corrections are however necessary to obtain correct result (|149|) . 

It seems plausible that we can add arbitrary function of the occupation numbers Si and to to T 
not spoiling the solution to the loop equation. However, if we assume normalizing conditions (|89|) 
and (|91j) . which follow from the condition of the "locality" of J-, that is, from the condition that J- 
depends only on local properties of dS near the branching points, we completely fix this ambiguity. 
The answer then becomes singular as Si — ► 0, and this is the singularity we are going to discuss now. 

The first type of singular behavior occurs when we shrink a cut with two "dynamical" branching 
points, e.g., bring fj,2 to //i. Setting ^ — f^i = e — > 0, we first consider what happens to the determinant 
term when closing the cycle A\. We have 



lim 



and 



hm <p — = <p ? * 1. 

10 This definitely implies the independence (vanishing) of _B-periods of dS under the action of d/dV(X). 
1 In |45) . the determinant term det<r is omitted from the answer. 
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Introducing 2/( n _i)(£) = — ^3) • • • (£ — M2n) , we then have for the determinant 

/ 1 X. £ £. £ 



lim cr, 



(n) 



V 7 ^ - ^3) • • • 0" - ^2n) 



det 



A 2 «-M)2/( TO -i)«) " ' JA 9 (€-Ai)»(„_l)(0 



n-2 



™ 2 «-*i)l/( n -i)(0 



™ 2 (£-/0»cn-i)(© 



*3 (?-/i)2/(n-l) (© 



A » (£-/0i/(n-l) (0 / 



and subtracting now the first column multiplied by <L 7 , 
= 2, . . . , g, we obtain just the determinant of (?( n -i) with the prefactor l/y( n _i) (/•*): 



from the fcth column for all 



lim det a 



M2 — > A i i=A i 



(n) 



det a 



(n-l)- 



(154) 



If we denote — [12 = e, then in the limit as /ii — > ^2 = 5*1 = e 2 2/(n-i)(A i )-^(/^) and the products 
of branching points are then combined exactly in a way to give 



lim Tf ] = F[ n ~ l) 



1 

12 



log5i + 0(e) 



(155) 



where we must substitute in the expression for ^ the new function M(£) = M (£)(£ — /-*) with the 
new zero at \i obtained by colliding the branching points [i\ and \i2- 

The singular part of T\ as Si — ► has therefore the structure 

^ n n— 1 

(^i)sing = -j^Y. l °Z Sl > S n = t -J2 s r ( 156 ) 

3=1 



i=i 



This is, in fact, a general phenomenon in the matrix model calculations: at any order of 't Hooft 
expansion, the only potential source for the singular contribution comes from degenerate geometry 
of curves, and it is related with the normalization factor in the matrix integral, that is, with the 
volume of (the orbit of) the unitary (sub)group. Indeed, the integral itself is a Taylor series in S^s 
(see formula (4.8) in w hhe the unitary group volume |51| in the case of Hermite polynomial 

distribution contributes the factor 



/ J] d\ k A 2 (A)e-*£ fc =i x l 
J -°° k=i 



Si/h s 2 
U r(Jfe) = e^- 

k=l 



G 2 (Si/h), 



(157) 



with G2(x) being the Barnes function |52j whose asymptotic expansion as x — ► 00 reads |52l I53j 

x_ 2 , 1 , 3,1.. ^ B* h 1 

2 



log G 2 (x) = — log x - — log x - -x 2 + -x log 2tt + C'(-l) + £ ^ _ ^ x2h _ 2 , 



(158) 
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where B 2 h are Bernoulli coefficients and £(s) is the Riemann (^-function. Discarding all polynomial 
terms, we obtain that the singular contribution as Si — * is |17| 

S 2 1 ^2h 1 

(^o) si ng= -y log 5i, (^l) 8 ing = -^ lo g^' Casing = ± h (h- 1) S 2h ~ 2 ' k ~ ^ ^ 

in the case where two dynamical branching points collide. This is the simplest way to segregate singular 
contribution; for instance, direct calculations demonstrate that the answer for genus-two correction in 
the one-cut case |2I.: has the prescribed singularity — 2iop-- 12 

In the hard-wall case, we also have the singularity appearing when a dynamical point collides with 
the hard wall, \x\ — > a%. In this case, we must calculate the asymptotic expressions for Laguerre, not 
Hermite, polynomials as the characteristic expression in this case is 

s 2 

r+00 Si / h t <! In / ft \ -4 Si 'h S 2 

/ J] dA fc A 2 (A)e-x^ = (Pi R J] (TO) 2 - e^ lo ^\G 2 (Si/h)f, (160) 
Jo k=i vtl/ fe=i 

because the norm for the normalized Laguerre polynomial -^r{x n e~ x ) is (n\) 2 . Then exactly the 
same considerations yield the corresponding singularities: 

1 2-£?2/i 1 

Casing = S i lo gSi, Casing = ~g lo § S *> (^)sing = 4^ _ fj ^-2 ' ^ > 2 - ( 161 ) 

Comparing with result (|149|) using the same considerations as above (in particular, the behavior of 
det a is exactly the same) , we find that 

lim =F[ n ~ X) -\\ogS x + 0{e), (162) 

where S\ ~ eM(//)i/( n _ 1 )(/i) and the function M(£) remains unchanged because no new zeros or poles 
arise when \i\ — ► a\. This demonstrates again the consistency of result (|149|) because the presence of 
additional factors depending on parameters ap in T\ would affect the (nonsingular) part of relation 

cm. 



5 Conclusion 

In conclusion, we briefly mention some related papers. Perhaps the first consideration of the WDVV 
equations for the differential dS with square-root singularities at the denominator was done from the 
matrix-model standpoint by Itoyama and Morozov |54| . There, it was used to regularize the matrix- 
model integral. Worth mentioning is also recent paper in which the CFT technique has been used 

12 In the proper normalization in which to = d 2 /16 for M(£) = 1. 
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for finding multiresolvents for the simplest, one-cut zero-potential matrix model with two hard edges. 
The asymptotic expansion for such an integral (the Legendre polynomials) reads 

log j a _ A 2 (A) n dX k = || log(o/8) + | log 2vr + 3C'(-1) - ~ log(t /h) 

k=l 

i « (h\ 2h - 2 B 2h ( 1 \ . , 

and in particular we see that the third derivative of the leading term vanishes, as follows from (|97l) . 
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